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Abstract 

We present a lattice calculation of the B s — B s transition matrix element 
through a four-quark operator Os = b(l — 75)5 6(1 — 75)5, which gives a 
leading contribution in the calculation of the width difference Ar s in the 
1/m.b expansion. The NRQCD formulation is used to describe b quark on the 
lattice. Using the next-to-leading formula of Beneke et at, we obtain (AT/r) s 
= 0.151(37)(45)(17), where the first error reflects the uncertainty of the B s 
meson decay constant, the second error comes from our calculation of the 
matrix element of O5, and the third represents unknown l/m& correction. 

PACS number(s): 12.38.Gc, 12.39.Hg, 13.20.He, 14.40.Nd 
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I. INTRODUCTION 



The mixing and decays of the B° s — B° s system play a complementary role to the B ± and 
B® — B% systems in studying the flavor mixing and CP violation [^j. In particular, if the 
width difference of the B° s — B° s system is sufficiently large, the angle ^3(7) of the unitarity 
triangle can be measured through untagged modes such as B s — > D^K^ or B s — > D*<p 
which would be promising not only because the method is theoretically clean but also 
feasible at future hadron colliders. 

The width difference AT b s of the B s — B s systems is calculated most reliably using 
the heavy quark expansion j|, and the size of a ratio (AT/T) Bs is roughly estimated as 
(Ar/r) Bs = 0.16(lo;o9). Now that the perturbative error has been reduced by the recent 
calculation of the next-to-leading order (NLO) QCD corrections ||, the largest remaining 
uncertainty comes from the matrix elements {B s \Ox{^b)\B s ) (X = L or S) of four-quark 
operators 

Ol = M 1 -75)sM 1 -T5)s, (1) 
O s = 6(l-75)s6(l-7 5 )s. (2) 

Lattice QCD is one of the most suitable tools for the nonperturbative computation of 
matrix elements such as the decay constants and the bag parameters. In fact a number 
of extensive studies, including ours 0, have already been done to obtain B^ j?J, which is 
a matrix element of the former operator Ol normalized by its vacuum saturation approx- 
imation. On the other hand, the matrix element Bs for the latter operator Os has been 
calculated in Ref. || only for the heavy-light meson around charm quark mass regime. It 
is required to perform a thorough study of Bs in order to give a reliable prediction of the 
B s width difference. The matrix element of Os is also required in the evaluation of the 
amplitude AM of the B^ — B^ mixing, if we assume the physics beyond the Standard 
Model such as the supersymmetric models |9]]. 

In this paper, we present a quenched lattice calculation of the matrix element of Os 
using the NRQCD formalism |TIJ for heavy quark and the 0(a)-improved Wilson action |TT 



for light quark. The NRQCD formalism is formulated as an inverse heavy quark mass 
expansion, and our action and operators consistently include entire 0(p/rriQ) terms, where 
p denotes a typical spatial momentum of a heavy quark inside a heavy-light meson. Higher 
order contribution of 0{p 2 /rnq) is also studied by introducing all necessary terms, and we 
find those effect is small for the b quark mass. 

In this work one-loop matching of the operator Os between continuum and lattice reg- 
ularizations is performed in the limit of infinitely heavy quark mass, so that the systematic 
error of 0(a s /(amQ)) is not removed. Since the b quark mass in the lattice unit is not 
extremely large, 0(a s /(arriQ)) gives a non-negligible effect in our final result, which could 
be as large as about 10% in a naive order counting argument. 

Using the NLO formula of Ref. || and our results for the matrix elements of Os as well as 
of Ol 0, we obtain a prediction (AT/T)b s = 0.151(37)(45)(17). The first error originates in 



1 We use the notation Bl instead of Bb to explicitly indicate that it represents a matrix element 
of O l . 
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the B s meson decay constant fs a = 245(30) MeV used to normalize the matrix elements, 
and the second is from our calculation of the matrix element of Os- The error from Bl is 
negligible, since it gives only a small contribution to the width difference. The last error is 
a crude estimate of the 0{l/m h ) correction as discussed in Refs. @]£J. 

This paper is organized as follows. We briefly summarize the NLO formula of Ref. || 
for the width difference in the next section. We present the perturbative matching of the 
operator Os in Section |TTTj, while the detail of the one-loop calculation is given in Appendix. 
We describe our simulation methods in Section [TV], and our results for the matrix element 
and the width difference are given in Section [V|. In section [VI|, we attempt to estimate the 
size of 0(a s / (am,Q)) error, which is specific to our work with NRQCD. Section [Vll| is devoted 
to a comparison of our result with a previous work by Gupta et al. |§, who obtained the 
same matrix element using the relativistic lattice action around charm quark mass. Finally 



our conclusion is given in Section |VIII| . A preliminary report of this work is included in 
Ref. WIl. 



II. WIDTH DIFFERENCE OF B s MESONS 

In this section we briefly summarize the formula to give the width difference of B s mesons, 
which was obtained by Beneke et al. in Ref. ||. 

The width difference in the B s — B s system is given by 

AT Ba = -2— l — (B s \Im i J d 4 xTH eff {x)H eff (0)\B s ), (3) 

where TC e ff is a AB=1 weak transition Hamiltonian. The main contribution comes from 
a transition bs — > cc followed by cc — > bs, and other contributions mediated by penguin 
operators are also included 0. 

Using the l/m b expansion, the transition operator Im i J d 4 xTTC e ff(x)Tt e ff(0) is repre- 
sented by the local four-quark operators Ol and 0g, which leads to the following formula 
at the next-to- leading order ||, 

(AT^ = 16n 2 B(B s ^Xev)f 2 B M B v 



V r J Bs g(z)r]Q C D m 



x [G{z)*~B L {m b ) + Gs(z)l^^ + VT^Tz5 1/m ^ . (4) 

Here, the quantity B(B S — > Xev) is the semi-leptonic decay branching ratio. The factors 
g(z) = 1 — 8z + 8z 3 — z 4 — 12z 2 lnz (z = ml/ml) and r/Q C D represent the phase space 
factor and the QCD correction, respectively. The coefficients G(z) and G${z) are functions 
including the next-to-leading QCD corrections, and their numerical values are given in Table 

1 of Ref. §. 

-Bf^mb) and Bs{rrib) are the B parameters defined with the MS scheme at the renor- 
malization scale fib = rrib- Their definitions are 

(B s \0 L (fi b )\B s ) (B s \0 L (fi b )\B s ) 



\(B s \A \0)(0\A \B s ) yiM\ 



d / \ _ (B a \O s (ji b )\B a ) _ (B s \O s (fib)\B s ) 

Bs ^ ] = »S»P^ - -IflMl X U ^ b) ' (6) 
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In the last expression in Eq. @, we change the normalization of {B s \Os{^ b )\B s ) with the 
decay constant fg a by factoring out the ratio 



(0\A o \B s ) 



(0\PQ*)\B.) 



(7) 



Using the equation of motion the ratio TZ(fi b ) is expressed in terms of the quark masses m b 
and m s as 



m b (fi b ) + m s (n h ) 



M, 



(8) 



where m b (/z 6 ) and m s (fj, b ) denote the quark masses defined with the MS scheme at scale \x b . 

Finally 8u m denotes l/m b corrections, which may be estimated using the factorization 
approximation [|]. 

Numerically evaluating the coefficients in the right hand side of Eq. (H), we obtain 



(AT 



B, 



245MeV 



0.008£ L (m ft ) + 0.150 



B s (m b ) 

n(m b y 



0.086 



(9) 



where we choose a recent world average of unquenched lattice simulations fs s = 245(30) 
MeV for the central value of the decay constant ]?]]. In the following sections we present a 
calculation of the parameter Bs(m b )/TZ(m b ) 2 . Our calculation of ^(mj,) is already available 
in Ref. HI. 



III. OPERATOR MATCHING 

In this section, we present the perturbative matching of continuum operator 0$ to the 
corresponding operators defined on the lattice. We follow the calculation method in Ref. 
|i~3|l , where the one- loop matching of the operator Ol is presented. 



Following the definition in Ref. [[J, we adopt modified minimal subtraction (MS) with 
the Naive Dimensional Regularization (NDR) scheme for the continuum operator Os{^ b ), 
in which 75 anticommutes with all The subtraction of evanescent operators is done 

with the definition given by Eqs. (13)-(15) of Ref. 0. The renormalization scale /i& is set to 
the b quark mass m b . 

While in the numerical simulations we apply the NRQCD formalism JTO] to the heavy 



quarks, in the perturbative calculation the heavy quarks are treated as a static quark [14 



More comments on this approximation will be given in the end of this section. The light 



quarks and gauge fields are described by the 0(a)-improved SW quark action |TTJ and the 
standard Wilson (plaquette) action, respectively, in both of the perturbative calculation and 
the numerical simulations. 

The operators involved in the calculation are 

O s = b l P L s l VP L s j , (10) 
6s = b l P LS i VP L s\ (11) 
O l = &VV Viols', (12) 
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Op = 2 I^Pls* V^PrJ + 4iV VPls* VP r s\ 
Or = &VV 

Old = ftVMT • «^)^ 

O pd = 2 6*7^(7 • aD)s i V^P L s j + AN b l P R tf ■ aD)s l VP L s j , 



(13) 
(14) 
(15) 
(16) 
(17) 



where Pl and Pr are chirality projection operators Pl/_r = 1^75- Color indices i and j run 
from one to N for SU(iV) gauge theory and a denotes lattice spacing. In the continuum, 
in which the chiral symmetry for light quark is preserved, the operator Os mixes only with 
Os and Ol- On the lattice, however, the chiral symmetry is explicitly broken with the 
SW action, so that additional operators with opposite chirality, Op and Or, appear in the 
operator matching. 

Other operators Osd, Old and Opd are higher dimensional operators introduced to 
cancel a discretization error of 0(a s a). However we neglect this discretization error in the 
numerical simulations. The result of 0(aa s ) matching coefficients is presented only for future 
use. 

Here we show the one-loop result of the matching. We leave the detail of the calculations 
for Appendix 0. The continuum operator Os(^b) is expressed by lattice operators O l x(l/a) 
as follows 



Os&b 



+s{H^+HS)-H]° 



lat i 
S 1 



I/O) 



a s 

47T 

47T 
47T 



3 

0.77 



Of(l/a) 



0.13 



Ojf ( 



-6.88 



2.58 



l/o) 



PE>(iA0 + 1 



1.15 



,.2,(1/0) 

O s 



(18) 



The operator Os is eliminated from the right hand side using a identity Os 
which is valid up to 0(p/mg). 

The heavy-light axial vector current A is also necessary to normalize the matrix element. 
The one-loop matching of A is already known as [p~6| , p~T| , p~3 



A = Z A {l/a)Af(l/a) + Z AD (l/a)A%(l/a 



Air 



2\n{a 2 m 2 b ) - 16.561 



K\l/a) 



4vr' 



13.01 ] Af Q (l/a) 



where A Q and A D0 are defined as 

A = &7o75S, 
-4do = &7o7s(7 -aD)s. 



(19) 



(20) 
(21) 



The higher dimensional operator Ap>o is introduced to remove the 0(a s a) errors. 

In Eqs. ( |i~8"D and (pT9|) , we apply the tadpole improvement [18|] using u — 1/8 A€ c clS clll 
average link variable. The normalization of the light quark filed is Jl — 3k/Ak c . 
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To obtain the matching coefficient for Bs/TZ 2 we combine Eqs. (|I8D and (|T9|), and linearize 
the perturbative expansion in a s . Omitting the higher dimensional operators, which we 
neglect in the following numerical simulations, we obtain 



B s (n b )/K(n b f 



4tt 1 a V 



a m 



+ 



Air 



In 



2 2 

a m b 




| +29.26 



5 



lat 



3.91 



5 



lat 



47T 



' 0.77 ~ 




0.13 









B 



lat 



where (X = 5, L, P, or i?) are '_B parameters' defined by 



B 



lat 
X 



(Of(l/a)) 



5 / Ala) 



(22) 



(23) 



(1/a)) 2 ' 

which we measure in the numerical simulations. 

Before closing this section, we should clarify the remaining uncertainty arising from the 
static approximation in the matching coefficients. In the simulation, the heavy quarks are 
described by the NRQCD action including the 0(p/mg) or 0(p 2 /mg) corrections consis- 
tently. The b quark field, which constitutes the operators measured in the simulation, is 
also improved through the same order as the action by the inverse Foldy-Wouthuysen-Tani 
transformation R^ 1 as 



R- 



Q 



where Q and x are the two-component quark and anti-quark fields in the NRQCD action. 
Therefore the truncation error only starts from 0(p 2 /rrvq) or 0(p 3 /wiq), which depends on 
the accuracy of our action and operators, even at the tree level matching. On the other hand, 
the static approximation in the perturbative calculation only leads to a lack of finite mass 
effects in the matching coefficients, but does not change the truncation error. Therefore, 
using the matching coefficients derived in this section the result has the 0(a s / (arriQ)) error. 



IV. SIMULATIONS 

The numerical simulations to extract B x ^ are almost the same as in our previous paper 
0, in which we calculated B L . We carried out a quenched simulation on 250 16 3 x 48 lattices 
at (3=5.9. The inverse lattice spacing from the string tension is 1.64 GeV. We employ the 
SW action for light quark with mean field improved c sw = 1/uq 3 with uq = 0.8734. The 



heavy quark is treated by two sets of NRQCD actions and fields |nj as was done in Ref . : 
one is truncated at 0{p/mQ) and the other includes entire 0(p 2 /tuq) corrections. We use 
the difference between the results from these sets to estimate the size of truncation error of 
the p/rriQ expansion. 

For the strong coupling constant used in the perturbative matching, we choose the V- 
scheme coupling ay(q*) with q* = 1/a, 2/a or it/ a. Their numerical values are cty(l/a) = 
0.270, a y (2/a)=0.193 and a v (ir/a) = 0.164. 

Other details of our simulations, such as the exact definition of the NRQCD action and 
the mass parameters used, are found in the previous paper 0. 
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V. RESULTS 



Figure [T] shows the mass dependence of B 1 ^ (X = S, L or P) defined in Eq. (|23f) . is 
equal to because of a symmetry under parity transformation. The light quark mass is 
interpolated to the strange quark mass. Since the light quark mass dependence is very small, 
in the following analysis we do not consider the error arising from the interpolation. The 
inverse heavy-light meson mass 1/Mp g , for which the light quark mass is also interpolated 
to the strange quark mass, is used as a horizontal axis. 

The difference between two results with different accuracies of the p/m,Q expansion does 
not exceed a few percent at the b quark mass, as explicitly presented in the figure by different 
symbols: circles for 0(p/rriQ) and triangles for 0(p 2 /itlq) accuracy. It justifies the use of 
the nonrelativistic expansion for the b quark. 

As we pointed out in the previous paper || , the vacuum saturation approximation ( VSA) 
gives a good approximation of the lattice data. In the static limit, it becomes Bg rSA ^=l, 
g(vsA) _ _g^/g an j B^SA) _ _g4^5_ p or finite heavy quark mass, the axial current and 
the pseudoscalar density involved in the VSA have different matrix elements. As a result, 
a mass dependence appears in the VSA of B x , as plotted by crosses (a flat line for B L ) 
in Fig. |l[ It is remarkable that the VSA explains the 1/Mp a dependence of the data very 
nicely. 

We combine the results for Bx to obtain Bs([Ji b ) /1Z(n b ) 2 using Eq. (|22|). The renormal- 
ization scale \i b is set to the b quark pole mass m b = 4.8 GeV according to Ref. ||. Figure || 
presents the 1/Mp s dependence of B s (rn b ) /1Z{m b ) 2 obtained with the 0(p/mg) (circles) and 
0(p 2 /mg) (triangles) accuracies and using ay(2/a)=0.193 as a coupling constant in the per- 
turbative matching. Typical size of the perturbative error may be evaluated by comparing 
the results obtained with different coupling constants. For this purpose, we also calculate 
the results with ay(7r/a)=0.164 and ay(l/a)=0.270, which are considered in the larger error 
bars in Fig. [|. We find that they give at most 5% differences at the b quark mass. 

Our numerical results interpolated to the physical B s meson mass Mb s = 5.37 GeV are 

• for 0{p/mo) accuracy 

\ 1.61(3) 



at 


q* 


= n/a 


at 


q* 


= 2/a 


at 


q* 


= 1/a 



(24) 



for 0{p 2 /rrig) accuracy 



B,(m h ) f L56 ( 3 ) at = 

-g^ = 1.59(3) at q* = 2/a , (25) 
K{m b ) [ L67 (3) at q * = i/ a 



where the error represents the statistical error. The variation due to the choice of the 
coupling constant av(q*) is explicitly shown. 

We attempt to estimate the size of systematic uncertainty in our result using an order 
counting of missing contributions. As we found in the previous paper |J, the dominant 
uncertainties are 
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0(a s /(am Q )) ~ 15%, 
0(a 2 s ) ~ 10%, 

0(a 2 Aj CJJ ) ~ 0(aA QC7D a s ) ~ 5%, 

when we assume Aq C d ~ 300 MeV and a s ~ 0.3. Although a naive order counting yields 
0(a s /(amQ)) ~ 10%, we take more conservative estimate ~ 15%, which is suggested in the 
study of bilinear operators as we will discuss in the next section. The effect of the truncation 
of the nonrelativistic expansion is negligible as we explicitly see in the difference between 
the two simulations of the 0(p/rriQ) and 0(p 2 /mn) accuracies. 
We finally obtain 

§g = 1.54(3)(30), (26) 

where the first error represents the statistical error, while the second is obtained by adding 
the sources of systematic uncertainty in quadrature. 

Using this result and the result for B L {m b ) previously obtained in Ref. [|], B L {m h ) = 
0.75(2)(12), we find 



= 0.151(37)(45)(17), (27) 

from Eq. The first error comes from the uncertainty in the decay constant fs s = 245(30) 
MeV, which is taken from the current world average of unquenched lattice calculations [|7|]. 
The second reflects the error in the calculation of Bs/TZ 2 presented above, and the last is 
obtained by assuming that the size of error in the l/m& correction Si/ m in Eq. (TJ) is ± 20%. 
The current experimental bound is (AT/T)b s < 0.42 |]T9 1 . 



The central value of our result in Eq. ( p7|) is much larger than the estimate 0.054^o;q32 
obtained by Beneke et al. ||. The main reasons are 

• The unquenched lattice result of fs s is about 15-20% larger than the previously known 
quenched result. 

• The central value of our result for Bs/TZ 2 is larger than the previous value obtained 
from the relativistic lattice calculation ||, which is used in Ref. ||. We will compare 



our result with theirs in section vTI. 



VI. FINITE MASS EFFECTS IN THE MATCHING COEFFICIENTS 

In this section, we attempt to estimate the size of 0(a s /(amQ)) error arising from the 
lack of necessary one-loop correction, by taking the ratio TZ(rrib) defined in Eq. (|7|) as an 
example. Although the 0(a s /(amg)) errors in bilinear operators and in the bag parameters 
are independent, it would still be useful to explicitly see the size of the error in a quantity, 
for which the correct one-loop coefficient is known. 

We compare the values of 7Z(rrib) 2 obtained with the following methods. 



S 



The quantity TZ(mb) 2 may be explicitly calculated in lattice simulation by measuring 
the matrix elements of axial-vector and pseudoscalar density. Results of the JLQCD 
collaboration obtained with the NRQCD action PO are plotted in Fig. ^ as a func- 



tion of 1/Mp s . One-loop matching to the continuum operator are calculated for two 
different lattice actions: static (filled circles) and NRQCD (open circles) f2T|j . 

2. The equation of motion may be used to obtain 

n( mb f = (™ b{mb) + ms{mb) ) 2 . (28) 



For the phenomenological values rhb(mb)=A.l-AA GeV and m s (2GeV)=0.06-0.17 
GeV P2"fl , which corresponds to m s (m fe )=0.05-0.14 GeV, we obtain 7?.(m&) 2 =0.66(5), 
which is shown by a star in Fig. |3|. 

The data obtained with the correct NRQCD matching coefficients (open circles) show a 
nice agreement with the phenomenological estimate (star). This suggests that the error 
in the calculation of the matrix element with correct matching coefficient is under good 
control. On the other hand, the data with the static matching coefficients (filled circles) are 
significantly lower, indicating a large systematic errors of 0(a s / '(amg)). The difference of 
TZ{mb) 2 between the two matching calculations is around 15% for the B s meson mass. We use 
this number for the estimation of the systematic error of 0[aj '{amq)) for Bs{rrib) /TZ(mb) 2 
in Sec. |V|. 



VII. DISCUSSION 

It is instructive to compare our result with the previous lattice calculation by Gupta, 
Bhattacharya and Sharpe ||, who used the Wilson fermion action for heavy quark with the 
mass around charm quark. Conversion of their result to the definition used in this paper is 
given in Ref. §, which yields 5 5 (2.33 GeV)=0.81 and 1^(2.33 GeV)=0.87. A 5-parameter 
for the operator Os (0) * s denoted as Bs- With the renormalization group evolution, it 
becomes Bs(rrib)=0.75 and Bs(rrib)=0.85 at fib = rrib- The error was not quoted except for 
the statistical one, which is 0.01 for each quantity. In order to compare the results obtained 
with different heavy quark mass, it is necessary to remove a logarithmic dependence on the 
heavy quark mass. We, therefore, define <&B s ( m b) as 



§B s {m h ) 



1-2 av "' In 



■2 



,a s {m b ) , frriQ\ n f N t 2a s (m b ) , ( m, 



5 Ait V mf 



B s(^b) + 7^ZT^ m ( Z% ) B s(m b ), (29) 



where rriQ denotes the heavy quark mass used in the simulation. In the calculation of Gupta 
et al. p| it is about the charm quark mass tuq = m c = 1.4 GeV. Using the coupling constant 
<y s { m b)= 0.22 corresponding to A^-=0.327 GeV and TZ{m b ) 2 obtained with the method 2 
in the previous section, we obtain 

®bT {m h )/K{m h ) 2 = 1.20, (30) 

which may be compared with our result of Bs(rrib) /TZ(rrib) 2 in Eq. (^6|). 
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The central value of our result is significantly higher than Eq. fl3~0|), which is one of 
the reasons of our larger value of compared to that of Ref. ||. We note, however, 

that the calculation with the unimproved relativistic action could suffer from large 0(arriQ) 
error, which is not even estimated in Ref. ||. In our NRQCD calculation, on the other hand, 
all possible systematic uncertainties are considered, but unfortunately the large systematic 
error of O (alpha s / amq) is left to be removed. Thus, at this stage we conclude that that the 
present accuracy of both calculations is not enough for a detailed comparison. 



VIII. CONCLUSION 

The width difference Ar s in the B s — B s mixing is expressed by the matrix elements 
of local four-quark operators in the l/m& expansion. The operator Os gives a dominant 
contribution among them and the nonperturbative calculation of its matrix element is es- 
sential for a reliable calculation of the width difference @||. We calculated a parameter 
Bsirrib) /7Z( m b) 2 , which is the matrix element normalized with a square of the B s meson 
decay constant as defined in Eq. (||), using lattice NRQCD formalism for heavy quark. 

From a quenched simulation at j3=5.9 with the 0(a)-improved light quark action, we 
obtain B s (mb)/lZ(mb) 2 = 1.54(3)(30), where statistical and systematic errors are given in 
that order. By explicitly performing two calculations with the different accuracies, we found 
that the 0(p 2 /ttlq) corrections in the NRQCD action and operators is only a few percent. 
One of the dominant sources of the systematic error is a lack of one-loop matching coefficients 
with finite mass corrections. We used the one-loop coefficients for the static action instead, 
which introduces a systematic error of order a s /(amo) ~ 15%. 

The large remaining uncertainty in our final result for (Ar/r) s , Eq. (P?D, comes partly 



from the error in our calculation of B 3(111^) / 1Z(mt) 2 . Another important source is present 
in the B s meson decay constant fs s , as it appears as j\ in the formula. 

We also discussed a comparison of our result with the previous one. We found that the 
central value of our result is significantly larger. However, since both calculations suffer from 
large systematic uncertainties, it would be fair to say that the discrepancy between the two 
results is not significant at the present level. 
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APPENDIX A: 

A matrix element of the continuum operator Os with free quark external states is ex- 
pressed at one-loop order as 
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1 + 



An 



13iV 2 - 18N + 9 -3iV 2 + 2N + 5 , 
1 — In 



AN 



2N 
3N 2 - AN 



+ 



a s 
An 



UN — 9 iV + 1 



2N 



N 



In 



+ 



mi 



A 2 



mi 



N 

2(N 



N 



\ mi 



(O s )o 
(O s )o 



a s 3(N-V 



+ 



An 

An 



AN 
2tt 1 

3 aX 



(Ol)o 
(@pd)o, 



(Al) 



where {Ox)o denotes a tree level matrix element of operator Ox, and the gluon mass A is 
introduced to regularize the infrared divergence. The evanescent operators are subtracted 
according to Eqs.(13)-(15) of Ref. ||. The expression is expanded in l/m b and only the 
leading terms are written. 



The corresponding expression for the lattice operator is | 23| ) |24 | 



(0f(l/a)> 



a s 
An 



-3N 2 + 2N + 5 



2N 
N 2 — 2 
N~ 



ln (a 2 A 2 ) + 



iV 2 -l 
2N 



(2) 



1 2N - 

d] c H 

2N 6N 



(v + v 1 )- 



N 



3N 



-Ji 



An 

a s 1 
H — - - 
An A 



N + l, ( 2 , 2 \ , 1 N-2, j N + 1 

In (a 2 A 2 ) + di + -c + ^ Ar (f + i/) + oAr Ji 



6iV 



3N 



(O s )o 
(O s )o 



d 2 - d 1 ] (Op) 



+ 
4 



a, N - 1 r 



4tt 2iV 
a s (iV + l)(iV 



u> + w 
2) 



+ 



47T 

An 2 



N 



(O R ) 

-(1 - c sw ) ln (a 2 A 2 ) - (V + F 7 )] (O SD )c 



a, 

47T 



(1 - c sw ) In (a 2 A 2 ) + (V + V 7 )] (O ld ) c 

T5T i O^O] <"»>.. 



(A2) 



where the constants c, d x , cf 2 , e > /, v i w d 1 , f 1 , v 1 , w 7 {7, f/ 7 , V, V 1 and J x are defined 
in Ref. [I4"|, 25| , p5 ||24",15| and their numerical values are tabulated in Table |. The coefficients 
with the superscript I denotes the terms appearing with the 0(a)-improvement. u comes 
from the tadpole improvement of the light quark wave function renormalization, and is also 
given in Table |I[ 

Matching above results and using a Fierz relation (Os)o = —(Os)o — 2 (Cl)o, which is 
satisfied in the static limit, we obtain for iV=3 



1 + S { 10 + l ln(aW ) + f lll (g)-i(/ + /' + e- + 1I f 
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— di + -c (v + v 1 ) + -Ji 

3 3 9 V 7 9 



+ - 



_«£ 1 

47T 4 

_«£ 1 r 
4vr 3 

Oi s 4 
+ — - 
4vr 3 

a s —1 
+ 4^ IT 



" 3 


--lnf 


2 2\ 






2 


3 V 









)f(l/a) 



d 2 -d'lo^l/a) 



u; + w 



0^(l/a) 



(1 - c sw ) In (a 2 A 2 ) + V + V 1 ] Of D (l/a) 



(1 - c sw ) In (a 2 A 2 ) + V + VjO^l/a) 



+- 



a, 1 



t/ + f/ 7 10^(l/a). 



47T 4 

A result with c sw =l is used in Eq. 



(A3) 
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TABLES 



c 


4.53 


d, 


5.46 


d 2 


-7.22 


d 1 


-4.13 


e (R) 


4.53 


f 


13.35 


f 1 


-3.64 


V 


-6.92 


v 1 


-6.72 


w 


-1.20 


w 1 


0.82 


Jl 


-4.85 


U 


4.89 


U 1 


-0.29 


V 


-7.14 


V 1 


1.98 


4 2) (lmk) 


-7T 2 


u ( q\k c ) 


-8.00 



TABLE I. Numerical values of parameters appearing in the one-loop lattice integrals. 
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FIGURES 
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FIG. 1. 1/Mp, dependence of (a) %\ (b) fijf (=5^*), and (c) B l p\ The results with 0(p/m Q ) 
accuracy (circles) are compared to those with 0(p 2 /itiq) (triangles) accuracy. The vacuum satu- 
ration approximation is shown by crosses. 
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FIG. 2. l/Mp s dependence of Bs(mb)/TZ(mb) 2 . Results with the 0{p/m,Q) (circles) and 
0{j)l thq) (triangles) accuracies are shown. The smaller error bars represent statistical errors, 
while the uncertainties, obtained from a quadratic sum of the statistical uncertainty and difference 
between the central values with ay(l/a) and ay(7r/a), are shown by the larger error bars. The 
central values are obtained with ay (2/ a). 




1/M p [GeV ] 



FIG. 3. 1/Mp s dependence of lZ(mb) 2 evaluated with the method 1 (circles) and 2 (star). See 
the text for the detail. Open and filled symbols are obtained with and without the l/(amg) 
corrections in the one-loop coefficients. 
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